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Introduction 


Lili  Chen,  Ph.D. 


This  project  is  aimed  at  exploring  MR  imaging  based  treatment  planning  for  radiotherapy  of  prostate 
cancer.  We  have  proposed  to  work  on  the  third  task  for  the  third  research  year.  The  tasks  include  (1) 
develop  practical  procedures  for  clinical  implementation  of  MRI  simulation;  (2)  develop  guidelines  for 
MRI-based  treatment  planning  dose  calculation,  and  (3)  develop  quality  assurance  programs  for  MRI 
simulation  for  prostate  cancer  treatment.  We  describe  our  work  for  the  third  year  as  follows. 

Body 

In  this  annual  report  we  report  on  the  research  accomplishments  associated  with  the  tasks  outlined  in 
the  approved  “Statement  of  Work”  task  3  between  Mar.  1,  2006  and  Feb.  28,  2007.  We  will  provide 
detailed  information  below  for  the  results  in  the  third  year. 

Task  3.  Develop  practical  procedures  for  clinical  implementation  of  MRI  simulation 

Creation  of  MR-based  Digitally  Reconstructed  Radiographs  (DRRs) 

During  this  period,  we  have  focused  on  developing  a  technique  to  create  MR-based  digitally 
reconstructed  radiographs  (DRR)  for  prostate  IMRT  patient  setup  when  MR-based  treatment  planning 
is  applied  clinically.  A  paper  entitled  “MRI-Based  Treatment  Planning  for  Prostate  IMRT:  Creation  of 
Digitally  Reconstructed  Radiographs  (DRR)”  has  been  accepted  by  Int.  J.  Radiat.  Oncol.  Biol.  Phys 
(Chen  et  al  2007a).  We  studied  MR  image  distortion  corrections  to  further  improve  the  accuracy  of 
dose  calculation  for  MR  based  treatment  planning  for  prostate  cancer.  A  paper  entitled  “Investigation 
of  MR  image  distortion  for  radiotherapy  treatment  planning  of  prostate  cancer”  has  been  published  in 
Physics  in  Medicine  and  Biology  (Chen  et  al  2006).  We  also  used  the  Monte  Carlo  method  to  further 
verify  dosimetric  accuracy  and  consistency  for  MR  based  IMRT  treatment  planning  for  prostate 
cancer.  A  short  paper  entitled  “Monte  Carlo  dose  verification  of  MR  image  based  IMRT  treatment 
planning  for  prostate  cancer”  has  been  submitted  to  XVth  Intenational  Conference  on  the  Use  of 
Computers  in  Radiation  Therapy  (ICCR.  Toronto,  Jun  2007)  (Chen  et  al  2007b).  The  two  manuscripts 
together  with  the  published  paper  are  also  attached  to  this  report.  We  summarize  the  results  and 
conclusions  of  these  studies  as  follows. 

Our  previous  studies  demonstrated  that  MRI-based  treatment  planning  meets  the  dosimetric  accuracy 
for  prostate  IMRT  and  it  is  adequate  to  use  unity  density  in  treatment  planning  dose  calculation  with 
co-planar  beam  arrangements  for  prostate  cancer  treatment  after  correction  of  MRI  distortions  (Chen  et 
al  2004a,  2004b).  With  CT-based  treatment  planning,  the  CT-based  DRRs  are  routinely  used  for 
patient  treatment  set-up  verification  by  comparing  with  portal  film  or  electronic  portal  imaging  devices 
(EPID).  With  MR-based  treatment  planning,  Since  MRI-derived  DRRs  do  not  provide  enough  bony 
structure  information  and  therefore  cannot  be  used  directly  for  checking  patient  positions.  To 
overcome  this  problem,  we  have  developed  a  technique  to  create  MR-based  DRRs  for  patient  initial 
setup  for  routine  clinical  applications  of  MR-based  treatment  planning  for  prostate  patient  treated  with 
IMRT.  Twenty  prostate  patients’  CT  and  MR  images  were  used  for  the  study.  CT  and  MR  images 
were  fused.  The  pelvic  bony  structures  including  femoral  heads,  pubic  rami,  ischium  and  ischial 
tuberosity  that  are  relevant  for  routine  clinical  patient  setup  were  manually  contoured  on  axial  MR 
images  using  the  AcQsim  planning  system.  The  contoured  bony  structures  were  then  assigned  a  bulk 
density  of  2.0  g/cm3.  The  MRI  based  DRRs  were  generated.  The  accuracy  of  the  MR  based  DDRs  was 
quantitatively  evaluated  by  comparing  MR-based  DRRs  with  CT-based  DRRs  for  these  patients.  For 
each  patient  8  measuring  points  on  both  coronal  and  sagittal  DRRs  were  used  for  quantitative 
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evaluation.  Our  results  showed  that  the  maximum  difference  in  the  mean  values  of  these  measurement 
points  is  1.3  and  the  maximum  difference  in  absolute  positions  is  within  3  mm  for  the  20  patients 
investigated.  MR-based  DRRs  are  comparable  to  CT-based  DRRs  for  prostate  IMRT.  This  technique 
has  been  used,  in  combination  with  the  BAT/in-room  CT  daily  target  localization  technique,  for  the 
clinical  implementation  of  MRI-based  treatment  planning  for  prostate  IMRT  at  FCCC. 

Develop  guidelines  for  MRI-based  treatment  planning  dose  calculation 

We  investigated  the  effect  of  MRI  residual  distortion  after  gradient  distortion  correction  (GDC)  on 
IMRT  treatment  planning  and  dosimetry  accuracy.  The  residual  distortion  errors  are  less  than  1  cm  and 
will  have  negligible  clinical  impact  for  more  than  90%  of  the  prostate  patients  whose  lateral 
dimensions  are  <40  cm.  We  have  investigated  on  MR  image  distortion  correction  to  further  improve 
the  accuracy  of  dose  calculation  for  MR  based  treatment  planning  for  prostate  cancer.  Our  studies 
showed  that,  with  our  routine  clinical  3-dimensinal  fast  spin  echo  sequences  (3DFSE,  256  x  256,  1.855 
mm  pixel,  TR  =  140  ms,  TE  =  3000  ms,  BW  readout  gradient  >  100  Hz/pixel),  there  was  no  patient- 
induced  susceptibility  distortions.  Therefore,  the  residual  machine  specific  geometrical  distortions  after 
the  GDC  could  be  quantified  by  phantom  measurements  and  further  reduced  by  our  point-by-point 
correction  technique.  The  effective  field  of  view  (FOVeff)  of  the  scanner  was  established  based  on  the 
actual  viewable  areas  with  adequate  geometric  distortion  corrections  (ensuring  <  5  mm  distortion 
error).  The  effective  FOVeff  for  prostate  imaging  using  a  standard  FOV  of  48  cm  has  been  expanded 
from  36  cm  using  the  existing  GDC  software  to  42  cm  using  the  point-by-point  distortion  correction 
technique  developed  in  this  work.  Our  results  indicated  that,  with  the  distortion  maps  established  in  this 
work,  we  could  correct  MR  geometrical  distortions  for  patients  of  lateral  dimensions  up  to  42  cm. 
Significant  improvement  in  dose  calculation  has  been  achieved  based  on  a  1-2  cm  improvement  in 
patient  external  contour  determination  (Chen  et  al  2006). 

Our  previous  study  results  showed  that  no  clinically  significant  differences  in  dose  calculations  were 
found  between  MRI-  and  CT-based  treatment  plans  using  the  same  beam  arrangements,  dose 
constraints  and  optimization  parameters.  We  also  validated  the  dosimetry  accuracy  of  MRI-based 
treatment  planning  by  recomputing  MRI-based  IMRT  plans  using  patient  CT  data  and  an  IMRT  QA 
phantom.  The  differences  in  dose  distributions  between  MRI  plans  and  the  corresponding  recomputed 
plans  were  generally  within  3%/3mm.  The  differences  in  isocenter  doses  between  MRI  dose 
calculation  and  phantom  measurements  were  within  our  clinical  criterion  of  4%. 

This  year  we  focused  on  the  Monte  Carlo  method  to  further  verify  dosimetric  accuracy  and  consistency 
for  MR  based  IMRT  treatment  planning  for  prostate  cancer.  The  Monte  Carlo  code  used  in  this  work 
was  MCSIM,  which  is  an  EGS4/PRESTA  user  code  developed  at  FCCC  (Ma  et  al  2002).  We  have 
performed  CT-based  IMRT  Monte  Carlo  dose  calculations  with  and  without  heterogeneity  corrections 
in  order  to  investigate  the  heterogeneity  effect  caused  by  different  beam  angle  arrangements.  Based  on 
the  results,  MR-based  IMRT  dose  calculations  were  performed  using  either  uniform  density  geometry 
or  unifonn  density  geometry  with  bulk  electron  density  assigned  to  bony  structures.  For  the  plans  with 
insignificant  inhomogeneity  effect,  uniform  geometries  with  water  density  were  used  in  the  MR-based 
dose  calculation.  For  the  plans  that  bony  structure  constitutes  a  large  part  of  volume  irradiated,  uniform 
density  geometry  with  bulk  electron  density  assigned  to  bony  structures  was  used  in  the  MR-based 
dose  calculation.  Each  IMRT  plan  was  evaluated  based  on  isodose  distributions  and  dose  volume 
histograms  (DVHs)  with  CT-based  or  MR-based  dose  calculations.  The  clinical  target  volume  (CTV) 
was  chosen  for  the  dose-volume  comparison.  Clinical  quantities  such  as  the  mean  dose,  maximum  and 
minimum  dose  received  by  the  CTV  and  the  critical  structures  were  compared.  The  maximum  dose 
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was  defined  as  the  highest  dose  received  by  1%  of  the  target  volume  and  the  minimum  dose  was 
defined  as  the  lowest  dose  received  by  99%  of  the  target  volume,  respectively.  Other  parameters  such 
as  the  dose  at  the  isocenter  and  the  dose  received  by  95%  and  5%  of  the  CTV  were  also  compared.  The 
paired  CT  and  MR  data  for  any  patients  in  this  work  were  pre-processed  to  have  the  same  pixel 
resolution.  The  internal  contours  of  the  targets  and  critical  organs  were  contoured  by  oncologists  on  the 
fused  CT-MR  images.  A  special  computer  code  was  developed  to  convert  the  patient  CT  and  MR 
image  data  from  the  DICOM  format  to  geometries  specially  formatted  for  the  MCSIM  code  (Chen  et  al 
2007). 

Our  results  showed  that  the  differences  in  dose  calculations  between  CT  data  (with  heterogeneity 
correction)  and  MR  data  (with  uniform  water  equivalent  geometry)  were  about  3%  or  less  and  less  than 
2%  in  the  mean  values  for  the  10  plans  with  beams  arranged  in  the  axial  plane.  For  MR-based 
calculations  (with  homogeneous  geometry),  our  results  demonstrated  that  the  differences  between  MR- 
based  calculations  and  CT-based  calculations  (without  heterogeneity  correction)  were  less  than  2%  for 
the  individual  patients  and  about  1%  in  the  mean  values,  which  proved  that  MR-based  IMRT  plans  can 
be  used  to  replace  CT-based  planning  clinically.  The  1%  -  2%  differences  in  dose  calculations  were 
mainly  caused  by  the  setup  uncertainties  of  the  two  imaging  modalities  if  geometrical  distortions  on 
the  MR  images  were  corrected  to  less  than  3  mm.  For  treatments  in  which  relatively  large  amount  of 
bones  are  irradiated,  MR-based  treatment  planning  with  homogeneous  geometry  would  not  be 
appropriate  because  of  the  excessive  attenuation  of  the  photon  beams  passing  through  bony  structures. 
However,  by  assigning  bulk  densities  to  the  bony  structures  especially  for  the  femurs  and  femoral 
heads,  the  dose  differences  could  be  reduced  to  less  than  3%.  The  bulk  density  assigned  to  the  femurs 
that  gave  the  best  fits  to  CT-based  calculations  with  heterogeneity  correction  was  1.8g/cm  in  our 
simulations  (Chen  et  al  2007). 

Develop  quality  assurance  programs  for  MRI  simulation  for  prostate  cancer  treatment 

We  have  established  a  practical  procedure  for  MR-based  treatment  planning.  1)  An  optimal  MR 
protocol  was  first  developed  for  contouring  the  target  and  critical  structures.  2)  We  investigated  the 
effect  of  MRI  residual  distortion  after  the  GDC  on  IMRT  treatment  planning  and  dosimetry  accuracy. 
The  residual  distortion  errors  are  less  than  1  cm  and  will  have  negligible  clinical  impact  for  more  than 
90%  of  the  prostate  patients  whose  lateral  dimensions  are  <40  cm.  For  patients  whose  lateral 
dimensions  are  >  40  cm  we  will  use  the  point-by-point  distortion  correction  technique  developed  in 
this  work.  Our  results  indicated  that,  with  the  distortion  maps  established  in  this  work,  we  could 
correct  MR  geometrical  distortions  for  patients  of  lateral  dimensions  up  to  42  cm.  3)  We  have 
investigated  optimal  fiducial  markers  for  MRI  simulation.  We  have  introduced  a  new  donut-shaped 
marker  to  improve  isocenter  definition  (IZI  medical  Product,  Baltimore,  Maryland  21244).  The  marker 
contains  iodine  with  a  1.5  cm  outer  diameter  and  a  4  mm  inner  diameter.  The  centers  of  the  markers 
can  be  detected  clearly  on  one  MR  slice  to  define  the  treatment  isocenter.  To  implement  MRI 
simulation,  a  set  of  trackable  lasers  has  been  installed  in  the  MR  room  for  patient  setup  and  isocenter 
determination.  4)  We  have  demonstrated  that  MRI-based  treatment  planning  meets  the  dosimetric 
accuracy  for  prostate  IMRT  and  it  is  adequate  to  use  unity  density  in  treatment  planning  dose 
calculation  with  co-planar  beams  for  prostate  cancer  treatment  after  MRI  distortion  corrections  (Chen 
et  al  2004a,  2004b),  5)  We  have  developed  practical  methods  for  heterogeneity  correction  for  MRI- 
based  dose  calculations  (Chen  et  al  2007).  6)  We  will  use  existing  Monte  Carlo  simulations  as  a  quality 
assurance  (QA)  tool,  implemented  at  FCCC,  for  IMRT  dose  verification  using  the  EGS4/MCSIM 
Monte  Carlo  program  (Ma  et  al  1999,  2000,  2002),  in  which  a  patient’s  CT  and  MLC  leaf  sequences 


6 


Lili  Chen,  Ph.D. 


are  directly  used  to  reconstruct  the  dose  to  be  received  by  the  patient.  We  will  MRI-scan  IMRT  QA 
phantoms  and  perforin  dose  calculations  to  determine  monitor  units  (MU)  for  IMRT  plans  and  then 
compare  with  ion  chamber  measurements;  7)  finally  MRI-based  DRRs  are  used  during  initial  treatment 
setup  together  with  CT-on-rails/cone-beam  CT//BAT  and  later  on  as  a  backup  for  these  imaging 
systems  if  the  systems  are  down. 


Key  Research  Accomplishments 

We  have  accomplished  the  following  tasks: 

•  The  dose  accuracy  for  MR-based  treatment  planning  of  prostate  cancer  has  been  validated 
using  the  Monte  Carlo  method  and  demonstrated  consistent  results. 

•  A  practical  method  for  heterogeneity  correction  for  MRI-based  dose  calculation  in 
inhomogeneous  patient  anatomy  has  been  developed. 

•  A  practical  technique  to  create  MR-based  DRRs  for  prostate  IMRT  has  been  developed  that  can 
be  used  for  patient  setup  when  MR-based  treatment  planning  is  applied  clinically. 


Reportable  Outcomes 

Peer-reviewed  papers  resulting  from  or  supported  in  part  by  this  grant : 

•  Chen  L,  Nguyen  T-B,  Jones  E,  Chen  Z,  Luo  W,  Wang  L,  Price  RA,  Pollack  A  and  Ma  C-M. 
MRI-Based  Treatment  Planning  for  Prostate  IMRT:  Creation  of  Digitally  Reconstructed 
Radiographs  (DRR).  International  Journal  of  Radiation  Oncology  Biology  Physics  (2007a) 
accepted 

•  Chen  Z,  Ma  C-M,  Paskalev  K,  Li  J,  Yang  J,  Richardson  T,  Palacio  L,  Xu  X  and  Chen  L. 
Investigation  of  MR  Image  Distortion  for  Radiotherapy  Treatment  Planning  of  Prostate  Cancer. 
Phys.  Med.  Biol.  51:1 393-1404  2006 

•  Chen  Z,  Ma  C-M,  Yang  J,  Li  J,  Luo  W,  Pan  J,  Paskalev  K.A,  Price  Jr  R  A,  Chen  Y  and  Chen 

L.  Monte  Carlo  dose  verification  of  MR  image  based  IMRT  treatment  planning  for  prostate 
cancer.  The  XVth  Internatinal  Conference  on  the  use  of  Computers  in  Radiation  Therapy, 
submitted  2007b 

•  Pan,  J.,  Li,  L.,  Chen,  L.,  Stathakis,  S.,  Luo,  W.,  du  Plessis,  P.,  Xiong,  W.,  Ma,  C.-M.  A 
Practical  Monte  Carlo  MU  Verification  Tool  for  IMRT  Quality  Assurance.  Phys.  Med.  Biol. 
51:2503-15,2006. 

•  Luo,  W.,  Li,  J.S.,  Price,  R.A.,  Chen,  L.,  Yang,  J.,  Pan,  J.,  Chen,  Z.,  McNeeley,  S.,  and  Ma,  C.- 

M.  Monte  Carlo  based  IMRT  QA  using  MLC  log  files  and  R/V  outputs,  Med.  Phys.  33:  2557- 
64,  2006. 
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•  Wang,  L.,  Li,  J.,  Paskalev,  K.,  Hoban,  P.,  Luo,  W.,  Chen,  L.,  McNeeley,  S.,  Price,  R.,  Ma,  C.- 
M.  Commissioning  and  quality  assurance  of  a  commercial  stereotactic  treatment  planning 
system  for  extracranial  IMRT.  JACMP,  7:  21-34,  2006. 

•  Wang,  L.,  Feigenberg,  S.,  Chen,  L.,  Paskalev,  K.,  Ma,  C.-M.  Benefit  of  3D  image-  guided 
stereotactic  localization  in  the  hypofractionated  treatment  of  lung  cancer.  Int.  J.  Rad.  Oncol. 
Bio.  Phys.  66:  738-747,  2006. 

Chapters: 

•  Y  Cao  and  Chen  L.  Current  MRI  Application  in  Radiation  Oncology.  AAPM  Summer  School 
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Conclusions 

We  have  made  significant  progress  during  our  third  year  investigation.  We  have  successfully 
performed  the  tasks  scheduled  in  the  “Statement  of  Work”.  We  have  developed  a  practical  technique  to 
create  MR-based  DRRs  for  prostate  IMRT  that  can  be  used  for  patient  setup  when  MR-based  treatment 
planning  is  applied  clinically.  We  have  developed  guidelines  for  MRI-based  treatment  planning  dose 
calculation  and  quality  assurance  programs  for  MRI  simulation  for  prostate  cancer  treatment. 
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Abstract 

Purpose:  Recent  studies  demonstrated  that  MRI-based  treatment  planning  meets  the 
dosimetric  accuracy  for  prostate  intensity-modulated  radiation  therapy  (IMRT)  and  it  is 
adequate  to  use  unity  density  in  treatment  planning  dose  calculation  for  prostate  after 
correction  of  MRI  distortions.  The  purpose  of  this  study  is  to  develop  a  technique  to  create 
MR-based  digitally  reconstructed  radiographs  (DRR)  for  patient  initial  setup  for  routine 
clinical  applications  of  MR-based  treatment  planning  for  prostate  IMRT. 

Materials  and  Methods:  Twenty  prostate  patients’  CT  and  MR  images  were  used  for  the 
study.  CT  and  MR  images  were  fused.  The  pelvic  bony  structures  including  femoral  heads, 
pubic  rami,  ischium  and  ischial  tuberosity  that  are  relevant  for  routine  clinical  patient  setup 
were  manually  contoured  on  axial  MR  images  using  the  AcQsim  planning  system.  The 
contoured  bony  structures  were  then  assigned  a  bulk  density  of  2.0  g/cm  .  The  MRI  based 
DRRs  were  generated.  The  accuracy  of  the  MR  based  DDRs  was  quantitatively  evaluated 
by  comparing  MR-based  DRRs  with  CT-based  DRRs  for  these  patients.  For  each  patient  8 
measuring  points  on  both  coronal  and  sagittal  DRRs  were  used  for  quantitative  evaluation. 
Results:  The  maximum  difference  in  the  mean  values  of  these  measurement  points  is  1.3 
±1.6  mm  and  the  maximum  difference  in  absolute  positions  is  within  3  mm  for  the  20 
patients  investigated. 

Conclusions:  MR-based  DRRs  are  comparable  to  CT-based  DRRs  for  prostate  IMRT  and 
can  be  used  for  patient  treatment  setup  when  MR-based  treatment  planning  is  applied 
clinically. 

Keywords:  Radiotherapy,  MRI  treatment  planning,  Prostate  Cancer,  MRI-based  DRRs 
Running  Title:  MRI-based  DRRs  for  Prostate  Cancer 


3 


Chen  4 


Introduction 

Magnetic  resonance  imaging  (MRI)  provides  superior  image  quality  for  soft  tissue 
delineation  over  computed  tomography  (CT)  and  is  widely  used  for  target  and  organ 
delineation  in  radiotherapy  for  treatment  planning  (1-4).  The  prostate  CTV  appears  to  be 
overestimated  with  CT  images  compared  with  MRI  (3,  5).  The  prostate  volume  on  CT 
appears  to  be  approximately  40%  larger  than  on  MRI  reported  by  Rasch  (5).  These  results 
were  consistent  with  those  reported  by  Krempien  et  al  (6).  Sannazzari  et  al  (3)  also  showed 
a  mean  overestimation  of  CTV  of  34%  with  CT  compared  with  MRI.  Results  showed  that 
radiation  doses  received  by  the  penile  bulb  and  corporal  bodies  may  cause  erectile 
dysfunction  after  prostate  treatment  (7).  Studies  are  being  performed  at  Fox  Chase  Cancer 
Center  (FCCC)  to  use  MRI  to  limit  the  dose  to  the  erectile  tissue  with  intensity-modulated 
radiation  therapy  (IMRT).  A  Phase  III  randomized  trial  at  FCCC  has  investigated  the 
clinical  significance  of  the  erectile  tissue-sparing  technique.  Recently,  Steenbakkers  et  al  (8) 
also  reported  that  the  dose  delivered  to  the  rectal  wall  and  bulb  of  the  penis  was 
significantly  reduced  in  treatment  plans  based  on  MRI-delineated  prostates  compared  with 
CT-delineated  prostates,  allowing  a  dose  escalation  of  2. 0-7.0  Gy  for  the  same  rectal  wall 
dose. 

Recent  studies  suggest  that  dose  escalation  with  3D  conformal  radiotherapy  (3DCRT)  or 
IMRT  can  potentially  improve  local  control  and  toxicity  (9-18).  Currently,  the  gold 
standard  for  prostate  treatment  planning  is  to  use  both  CT  and  MR  imaging  to  take  the 
advantage  of  the  electron  density  information  with  the  CT  and  the  superior  soft  tissue 
capabilities  with  MRI.  Ideally,  treatment  planning  would  be  based  solely  on  a  clinically 
useful  diagnostic  MR  study.  The  advantages  of  MR-based  treatment  planning  are  1)  to 
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remove  potential  errors  due  to  image  fusion;  2)  to  reduce  patient  cost  by  avoiding 
redundant  CT  scans;  and  3)  to  eliminate  unnecessary  radiation  exposure  to  the  patient  and 
to  save  patient,  staff  and  machine  time.  However,  there  are  several  challenges  to  directly 
using  MR  images  for  treatment  planning,  including  (1)  the  lack  of  electron  density 
information  that  is  needed  for  heterogeneity  corrections  in  dose  calculation,  (2)  image 
distortion  that  affects  patient  external  contour  determination  and  therefore  introduces  dose 
calculation  uncertainty,  and  (3)  lack  of  bony  structures  to  derive  effective  DRRs  for  patient 
setup.  The  first  two  issues,  i.e.  the  lack  of  electron  density  information  and  image 
distortions  in  MR-based  prostate  treatment  planning  have  been  discussed  in  great  detail  by 
Chen  et  al  (19).  The  purpose  of  this  study  is  to  investigate  a  technique  to  create  MR-based 
DRRs  for  prostate  IMRT  patient  setup  when  MR-based  treatment  planning  is  applied  clinically. 

Material  and  Methods 
Image  collection 

CT-MR  fusion  with  CT-based  dose  calculation  has  been  a  routine  procedure  for  prostate 
IMRT  treatment  planning  at  FCCC  since  2001.  CT-based  DRRs  together  with  BAT  (B- 
mode  Acquisition  and  Targeting,  NOMOS,  Corp.,  Sewickley,  PA)  ultrasound  or  the  CT- 
on-rails  system  (Siemens  Medical  Solutions,  Concord,  CA)  are  used  for  initial  patient  setup 
and  target  localization.  On-line  images,  obtained  from  either  the  CT-on-rails  system  or  the 
BAT  ultrasound  system  were  used  for  daily  target  localization  to  correct  for  prostate  inter¬ 
fraction  motion.  For  this  study,  twenty  patients  with  fused  CT-MR  images  were  used. 
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Calibration  of  bony  structures  on  MRI 

Figure  1  shows  an  example  of  the  DRR  directly  created  from  MR  images.  It  is  clearly 
demonstrated  that  the  DRRs  do  not  provide  enough  bony  structure  information  for  patient 
positioning  verification.  Therefore,  it  is  necessary  to  develop  a  technique  for  the  creation  of 
useful  MR-based  DRRs  for  initial  treatment  setup  when  CT  is  not  involved.  Clinical  bony 
land  marks  such  as  pubis,  ischium  and  ischial  tuberosity  on  CT-based  DRRs  are  used  to 
adjust  the  patient  treatment  position  by  comparing  with  anterior-posterior  (AP)  and  lateral 
portal  images.  Our  goal  is  to  develop  a  technique  to  create  MR-based  DRRs  that  are 
comparable  to  CT-based  DRRs  for  the  clinical  implementation  of  MR-based  treatment 
planning  for  prostate  IMRT.  Our  approach  is  to  generate  MR-based  DRRs  by  contouring 
the  relevant  bony  structures  and  assign  them  with  a  bulk  density.  We  perfonned  this  study 
on  an  AcQsim  system  (Philips  Medical  System,  Cleveland,  OH). 


On  T2  weighted  MR  images,  bones  appear  hyperintense  (white)  surrounded  by  narrow 
hypo-intense  (dark)  regions.  The  marrow  inside  the  bone  or  soft  bones  shows  high  intensity 
(white)  while  compact  bone  shows  lower  intensity  (dark).  Therefore,  the  boundary  between 
the  bone  and  the  soft  tissue  is  not  always  distinguishable  on  MRI.  In  order  to  study  the 
correlation  between  CT  number  (electron  density)  and  MR  signal  intensity,  CT  and  MR 
images  were  fused  according  to  bony  anatomy  using  the  interactive  method  available  on 
the  AcQsim  system.  CT  images  were  loaded  as  primary  images  and  MR  images  were 
loaded  as  secondary  images.  The  two  images  were  shown  side  by  side  on  the  same 
computer  screen.  The  bones  were  contoured  based  on  CT  images  while  the  contours  were 
also  shown  on  MRI  simultaneously.  The  calibration  of  the  bones  on  MRI  that  are  relevant 


6 


Chen  7 


to  prostate  treatment  setups  was  established  based  on  CT  images.  Based  on  this  calibration, 
the  boundary  of  each  relevant  bone  can  be  easily  determined  on  MRI  alone. 

Guidelines  for  bone  contouring 

The  technique  that  we  describe  here  is  intended  to  serve  as  a  guide  for  the  clinical 
implementation  of  MR-based  treatment  planning  for  prostate  cancer.  We  outline  contours 
of  the  bones  at  the  level  of  the  top  of  the  femoral  head  (a)  to  the  lesser  trochanter  of  the 
femur  (b)  as  shown  in  Figure  2.  In  order  to  describe  our  contour  technique  clearly  we 
categorize  the  pelvic  bones  into  three  groups  according  to  the  similarity  of  the  bone  shape 
as  shown  in  Figure  2.  Figure  3  shows  the  correlation  of  the  bony  contours  on  the  axial  MR 
image  and  contour  lines  on  the  coronal  and  sagittal  MR  images.  Critical  points  for  each 
group  are  determined  based  on  the  relation  of  the  points  on  the  DRR  that  are  used  for 
checking  patient  positions  routinely. 

For  convenience,  we  have  defined  a  few  technical  terms  for  our  contouring  technique 
including  “white  side  (strict)”,  “white  side  (non-strict)”,  “dark  side  (strict)”,  and  “dark  side 
(non-strict)”.  The  term  “white”  indicates  the  hyper-intense  region  on  the  T2  weighted 
image  and  “dark”  indicates  the  hypo-intense  region  on  the  T2  weighted  image.  Thus, 
“white  side  (strict)”  means  to  draw  a  contour  immediately  outside  the  hyper-intense  region 
and  “white  side  (non-strict)”  means  the  contour  is  drawn  slightly  away  from  the  hyper¬ 
intense  region  (about  1  pixel  away),  as  shown  in  Figures  4.  The  same  concept  can  be 
applied  to  the  terms  “dark  side  (strict)”  and  “dark  side  (non-strict)”.  As  mentioned 
previously  we  divided  the  pelvic  region  into  3  groups.  For  each  group  the  left  and  right 
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bones  are  mirror  images  on  the  axial  view.  Therefore  we  only  describe  our  method  for  the 
right-side  bones.  To  start  contouring,  “Virtual  fluoroscopy”  is  selected  on  the  tool  menu, 
which  displays  the  three  dimensional  views  of  the  MR  images.  The  contours  are  drawn  on 
the  axial  images  slice  by  slice  while  the  contours  will  show  simultaneously  on  both  coronal 
and  sagittal  views. 

Contouring  bony  structures  for  the  first  group  (region  1) 

The  points  of  interest  (POI)  here  are  “A”  and  “B”  (see  Figure  5)  that  will  determine  points 
“a”  and  “b”  on  the  coronal  DRR  and  “a”  on  the  sagittal  DRR  (see  Figure  9).  Point  “C”  was 
found  not  critical  for  patient  setup  that  will  appear  as  “b”  on  the  sagittal  DRR. 

To  contour  POI  A,  one  can  imagine  that  there  are  2  lines  “a”  and  “b”  on  two  sides  of  the 
peak.  The  intersection  of  “a”  and  “b”  is  the  POI  A.  Lines  “a”  and  “b”  are  determined  based 
on  the  boundary  between  hyper-  and  hypo-intense  regions  using  the  “white  side  (non- 
strict)”  technique.  POI  B  is  the  intersection  of  the  straight  lines  between  black  side  and 
white  side  as  shown  in  Figure  6.  It  should  be  mentioned  that  the  exact  locations  of  the  POI 
points  will  affect  the  accuracy  of  the  resulting  DRR  while  the  volumes  of  the  contoured 
structures  will  not. 

Contouring  bony  structures  for  the  second  group  (region  2) 

There  are  two  relevant  bones  in  region  2  including  pubic  symphysis  and  ischial  tuberosity. 
As  shown  in  Figure  7,  POI  A  and  POI  C  are  critical  for  the  positioning  check  clinically 
while  POI  B,  E  and  D  are  not  critical.  For  point  A  the  vertical  line  is  contoured  with  “white 
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side  (strict)”.  Attention  should  be  paid  to  the  horizontal  line  here  because  its  height  will 
appear  on  the  sagittal  DRR  that  will  be  used  to  judge  the  patient  AP  position.  The 
horizontal  line  is  contoured  with  “white  side  (strict)”.  We  draw  the  horizontal  line  starting 
from  where  the  white  (hyper-intense)  begins  to  blur.  For  POI  C  we  use  the  same  method  as 
for  POI  B  of  the  first  group.  For  POI  B,  E  and  D  we  contour  with  “white  side  (non-strict)” 
since  they  are  not  critical  to  patient  setups. 

Contouring  bony  structures  for  the  third  group  (region  3) 

For  the  third  group,  only  POI  A  is  the  critical  point.  The  way  to  determine  this  point  is 
similar  to  that  for  POI  A  for  the  first  group.  This  means  that  POI  A  is  the  intersection  of 
two  lines:  one  is  on  top  (white  side)  and  the  other  is  the  dark  side  (non-strict)  of  the  gap  as 
shown  in  Figure  8. 

Creation  of  MR-based  DRRs 

After  contouring  the  three  bone  groups,  proper  names  were  given  to  these  bones  in  the 
“organ  selection”  menu  (e.g.,  symphysis).  A  bulk  electron  density  was  assigned  to  these 
bones.  Based  on  our  experience  an  electron  density  between  1.8  and  2.0g/mm  is  adequate. 
After  assigning  electron  densities  to  all  the  relevant  bones  the  DRRs  were  generated.  The 
quality  of  the  coronal  and  lateral  DRRs  can  be  improved  by  adjusting  the  window  and  level 
parameters.  Figure  9  shows  the  MR-based  DRR  and  the  CT-based  DRR  for  one  of  the  20 
patients  investigated. 


9 


Chen  10 


Data  analysis 

Paired  MR-based  and  CT-based  DRRs  were  compared  quantitatively  for  all  the  patients 
investigated.  Contoured  bony  structures  from  MRI  were  superimposed  on  CT-based  DRRs 
using  the  Photoshop  software  with  the  same  isocenter  (Figure  10).  The  8  points  from  both 
coronal  and  sagittal  DRRs  were  used  for  the  comparison.  The  difference  for  each  position 
point  between  CT  and  MRI  was  calculated  based  on  the  pixel  value.  Minus  values  mean 
the  position  points  on  MR-based  DRRs  were  inside  the  outlines  of  the  corresponding  bony 
structures  on  the  CT-based  DRRs  while  positive  values  mean  the  MR-determined  position 
points  were  outside  the  outlines  of  bony  structures  on  the  CT-based  DRRs.  The  mean  and  1 
standard  deviation  of  the  error  for  each  measurement  point  were  calculated  for  20  patients. 
The  maximum  differences  between  CT  and  MR  for  each  measurement  point  were  also 
calculated. 

Results  and  discussions 

In  this  work,  we  have  developed  a  method  to  generate  MR-based  DRRs  for  the  clinical 
implementation  of  MRI-based  treatment  planning  for  prostate  IMRT.  Figure  9  shows  MR- 
based  DRRs  and  CT-based  DRRs  for  one  of  the  patients  investigated  in  this  work.  The 
MR-based  DRRs  were  generated  with  a  2.0  g/cm  density  assigned  to  the  relevant  bony 
structures.  It  demonstrates  that  MR-based  DRRs  are  comparable  to  CT-based  DRRs  and  it 
is  reasonable  to  use  2.0  g/cm3  as  the  bulk  density.  In  fact,  one  can  also  select  to  print  MR- 
based  DRRs  with  the  outlined  bony  contours  on  top  of  the  gray-scale  image.  This  way,  the 
contour  lines  (instead  of  the  gray-scale  image)  on  the  DRRs  can  be  directly  used  to 
compare  the  bony  structures  on  a  portal  film  or  an  EPID  image  (see  Figure  10). 
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To  evaluate  the  accuracy  of  the  MR-based  DRRs  generated  using  our  method,  we  have 
selected  8  relevant  position  points  on  coronal  and  sagittal  DRRs  to  compare  the  actual 
locations  of  the  relevant  bony  structures.  The  differences  of  the  position  points  between 
MR-based  DRRs  and  CT-based  DRRs  for  20  patients  are  summarized  in  Table  1.  Since  the 
patients’  CT  and  MR  images  were  directly  taken  from  our  clinical  data  base  and  they  were 
already  fused  for  treatment  planning  the  positional  differences  shown  here  also  included 
the  positional  errors  due  to  the  fusion  process,  which  is  estimated  to  be  about  2mm.  The 
agreement  between  MR-based  DRRs  and  CT-based  DRRs  is  excellent.  The  mean 
difference  is  within  0.8  ±1.5  mm  except  for  point  7  (1.3  ±  1.6)  mm  and  the  maximum 
difference  is  within  3.0  mm  for  all  180  position  points  evaluated  for  the  20  patients  except 
for  point  7  of  patient  19  (3.6mm).  From  a  clinical  point  of  view,  point  7  is  less  important 
than  point  8  (pubis)  that  is  used  for  the  determination  of  the  patient  AP  positions.  These 
results  demonstrate  that  there  is  no  systematic  biasing  in  our  contouring  method. 

It  should  be  noted  that  on  the  coronal  DRR  all  the  bones  (and  position  points)  have  mirror 
images  and  therefore  their  relative  positions  are  less  affected  by  the  contouring  method 
used  as  long  as  they  are  contoured  consistently.  For  example,  if  the  symphysis  is  contoured 
consistently  between  the  left  side  and  the  right  side  the  midline  of  the  symphysis  (i.e.,  the 
left-right  position  of  the  patient)  determined  by  points  3  and  4  will  be  accurate  even  though 
the  absolute  positions  of  points  3  and  4  may  be  affected  by  the  contouring  method  or  by 
inter-personal  variations.  The  same  can  be  said  about  the  patient  superior-inferior  position 
which  can  be  accurately  determined  by  bony  landmarks  such  as  pubis  and  ischium  through 
paired  points  such  as  (1,  5)  and  (2,  6).  The  patient  AP  position  is  mainly  determined  based 
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on  the  pubis  position  (i.e.,  point  8  on  the  sagittal  DRR).  Again,  the  relative  position 
determined  by  a  bony  structure  is  expected  to  be  more  accurate  than  the  absolute  positions 
of  the  contour  points  defining  the  bony  structure.  The  same  phenomenon  can  be  observed 
with  gray-scale  images  when  different  “window”  and  “level”  parameters  are  used.  Based 
on  this  analysis  we  estimate  the  accuracy  of  our  MR-based  DRRs  is  about  3  mm  or  better. 
This  does  not  include  the  uncertainties  introduced  by  the  CT-MRI  fusion  process  since  they 
will  be  eliminated  when  MRI  alone  is  used  for  radiotherapy  treatment  planning. 

MR-based  DRRs  can  be  used  for  initial  patient  position  verification  by  comparing  with 
images  taken  using  portal  film  or  an  EPID.  In  our  clinical  implementation  of  MRI-based 
treatment  planning  for  prostate  IMRT,  MRI-based  DRRs  are  used  during  initial  treatment 
setups  together  with  either  a  BAT  or  a  CT-on-rails  system  to  set  the  baseline  treatment 
position.  The  3  mm  accuracy  for  MR-based  DRRs  is  considered  adequate  for  this  purpose. 
In  fact,  for  routine  prostate  IMRT  treatment  at  FCCC,  we  use  the  BAT  ultrasound  system 
or  the  CT-on-rails  system  as  a  primary  means  to  relocate  the  soft  tissue  target  prior  to  each 
treatment.  The  initial  patient  setup  based  on  portal  film  or  EPID  only  serves  as  a  secondary 
means  to  place  the  treatment  target  0.5  -  0.7  cm  within  the  treatment  isocenter.  This  will 
facilitate  the  target  localization  process  in  the  subsequent  daily  treatments.  Table  2  shows 
the  mean  shift  values  of  the  38  treatment  fractions  for  the  first  4  prostate  IMRT  patients 
planned  using  MRI  alone  and  4  other  prostate  IMRT  patients  planned  using  CT  alone, 
being  treated  at  the  same  time.  The  shift  values  in  the  three  major  axes  were  all  within 
expected  values.  Occasionally  larger  shifts  (>lcm)  were  observed  for  both  CT-  and  MRI- 
based  planning  when  large  gas  pockets  occurred  in  the  rectum.  The  combined  target 
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localization  accuracy  for  routine  prostate  treatment  is  expected  to  be  better  than  5  mm  at 
95%  confidence  level. 

Finally,  we  also  investigated  the  time  required  to  generate  MR-based  DRRs.  Although 
manual  contouring  is  needed  for  making  practical  MR-based  DRRs  the  time  required  to 
contour  all  relevant  bones  for  a  patient  was  approximately  10-15  minutes.  The  tune  to 
assign  a  bulk  electron  density  to  the  relevant  bony  structures  was  minimal.  In  fact,  it  is  not 
necessary  to  contour  the  femoral  head  since  the  pubis  (symphysis),  ischium  and  ischial 
tuberosity  positions  are  sufficient  to  determine  the  patient  positions.  We  have  contoured 
femoral  heads  to  provide  redundant  checks  and  more  realistic  DRRs. 

In  summary,  we  have  explored  the  use  of  MR-based  DRRs  for  initial  radiotherapy 
treatment  setup  verification  for  prostate  IMRT  by  comparing  with  images  taken  using 
portal  film  or  EPID.  Our  results  showed  that  MR-based  DRRs  utilizing  the  outlines  of 
relevant  bony  structures  have  an  accuracy  of  about  3  mm,  which  are  adequate  for  initial 
patient  setup.  This  technique  has  been  used,  in  combination  with  the  BAT/in-room  CT 
daily  target  localization  technique,  for  the  clinical  implementation  of  MRI-based  treatment 
planning  for  prostate  IMRT  at  FCCC  since  Nov.  2003. 
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Figure  Legends 

Figure  1.  Digitally  reconstructed  radiographs  (DRR)  from  MRI  (a)  and  CT  (b).  Top  figures 
are  coronal  view  and  bottom  figures  are  sagittal  view  of  the  pelvic  region. 

Figure  2.  Delineation  of  bony  structures  for  MR-based  DRRs  showing  contours  starting  at 
the  level  of  the  femoral  head  (a)  to  the  lesser  trochanter  of  the  femur  (b).  The  pelvic 
bones  are  divided  into  three  groups  according  to  the  similarity  of  the  bone  shapes. 

Figure  3.  Correlations  of  the  relevant  bone  contours  on  axial  and  coronal  views  (a)  and 
axial  and  sagittal  views  (b).  Detailed  descriptions  are  given  in  the  text. 

Figure  4.  Demonstration  of  the  contouring  technique:  (a)  “Contour  White  side  (strict)” 
means  contouring  immediately  outside  the  hyper-intense  region,  (b)  “Contour  white 
side  (non-strict)”  means  contouring  slightly  away  from  the  hyper-intense  region  (by 
about  1  pixel). 

Figure  5.  An  example  of  the  first  bone  group.  POI  A  is  correlated  to  both  the  coronal  and 
sagittal  DRRs  while  POI  B  is  correlated  to  the  coronal  DRR  and  POI  C  is  correlated 
to  the  sagittal  DRR  only. 

Figure  6.  An  example  to  delineate  POI  A  and  POI  B:  they  are  the  intersection  points  of  the 
lines  “a”  and  “b”  (a);  The  volume  defined  by  the  contours  away  from  the  POIs  will 
not  affect  the  accuracy  of  the  DRR  (b). 

Figure  7.  An  example  of  the  second  bone  group  and  the  correlation  of  the  axial  contours 
with  the  coronal  and  sagittal  DRRs. 

Figure  8.  An  example  of  the  third  bone  group  and  the  correlation  of  the  axial  contours  with 
the  coronal  and  sagittal  DRRs. 
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Figure  9.  An  example  of  MR-based  DRRs  (a)  and  CT-based  DRRs  (b).  The  MR-based 
DRRs  were  generated  after  the  relevant  bones  have  been  contoured  and  assigned  to 
a  bulk  density  of  2.0  g/cm3. 

Figure  10.  Measurement  positions  on  coronal  and  sagittal  DRRs  to  evaluate  the  accuracy  of 
MR-based  DRRs  (characterized  by  the  MR-based  outlines/contour  lines).  Points  3 
and  4  are  used  for  the  L-R  position  setup,  points  1  and  5  or  2  and  6  are  used  for  the 
I-S  position  setup,  and  points  7  and  8  are  used  for  the  A-P  position  setup. 
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Table  2  Mean  shift  values  along  the  three  major  axes  based  on  the  BAT  ultrasound  target 
localization  procedures  for  38  IMRT  treatments.  The  shift  values  were  calculated  as  the  differences 
between  the  planned  target  locations  and  the  actual  target  locations  measured  by  the  BAT  system 
prior  to  each  treatment.  The  initial  treatment  setups  were  using  either  CT-based  DRRs  or  MR-based 
DRRs. 


Patient  # 

Shifts  based  on  MR-simulation  (cm) 

RL 

AP 

SI 

1 

-0.19  ±0.25 

0.07  ±0.39 

0.51  ±0.43 

2 

-0.20  ±0.32 

0.31  ±0.53 

0.10  ±0.48 

3 

0.15  ±0.25 

0.39  ±0.23 

0.78  ±0.23 

4 

0.12  ±0.29 

0.53  ±0.29 

0.14  ±0.29 

Patient  # 

Shifts  based  on  CT- simulation  (cm) 

RL 

AP 

SI 

1 

0.09  ±0.42 

-0.33  ±0.30 

-0.03  ±0.34 

2 

0.28  ±0.25 

0.03  ±  0.42 

0.48  ±0.41 

3 

-0.01  ±0.25 

-0.06  ±0.27 

0.32  ±0.32 

4 

0.05  ±0.24 

0.04  ±0.17 

0.40  ±0.40 
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MR  (a)  CT  (b) 


Figure  1.  Digitally  reconstructed  radiographs  (DRR)  from  MRI  (a)  and  CT  (b).  Top 
figures  are  coronal  view  and  bottom  figures  are  sagittal  view  of  the  pelvic  region. 
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Figure  2.  Delineation  of  bony  structures  for  MR-based  DRRs  showing  contours 
starting  at  the  level  of  the  femoral  head  (a)  to  the  lesser  trochanter  of  the  femur  (b). 
The  pelvic  bones  are  divided  into  three  groups  according  to  the  similarity  of  the 
bone  shapes. 
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(a) 


(b) 


Figure  3.  Correlations  of  the  relevant  bone  contours  on  axial  and  coronal  views  (a)  and  axial  and 
sagittal  views  (b).  Detailed  descriptions  are  given  in  the  text. 
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(b) 

Figure  4.  Demonstration  of  the  contouring  technique:  (a)  “Contour  White  side  (strict)”  means 
contouring  immediately  outside  the  hyper-intense  region,  (b)  “Contour  white  side  (non-strict)” 
means  contouring  slightly  away  from  the  hyper-intense  region  (by  about  1  pixel). 
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Figure  5.  An  example  of  the  first  bone  group.  POI  A  is  correlated  to  both  the  coronal 
and  sagittal  DRRs  while  POI  B  is  correlated  to  the  coronal  DRR  and  POI  C  is 
correlated  to  the  sagittal  DRR  only. 
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Figure  6.  An  example  to  delineate  POI  A  and  POI  B:  they  are  the  intersection  points  of 
the  lines  “a”  and  “b”  (a);  The  volume  defined  by  the  contours  away  from  the  POIs  will 
not  affect  the  accuracy  of  the  DRR  (b). 

rigure  a.  /\ri  example  or  me  msi  Done  group,  r kj i  /\  is  coireraieu  lo  oom  me  eoronar 
and  sagittal  DRRs  while  POI  B  is  correlated  to  the  coronal  DRR  and  POI  C  is 
correlated  to  the  sagittal  DRR  only. 
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Contour  the  second  group 


◄ - coronal  view 


Figure  7.  An  example  of  the  second  bone  group  and  the  correlation  of  the 
axial  contours  with  the  coronal  and  sagittal  DRRs. 
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Figure  8.  An  example  of  the  third  bone  group  and  the  correlation  of  the  axial  contours 
with  the  coronal  and  sagittal  DRRs. 
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Figure  9.  An  example  of  MR-based  DRRs  (a)  and  CT-based  DRRs  (b).  The  MR- 
based  DRRs  were  generated  after  the  relevant  bones  have  been  contoured  and 

•5 

assigned  to  a  bulk  density  of  2.0  g/cnr . 
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Figure  10.  Measurement  positions  on  coronal  and  sagittal  DRRs  to  evaluate  the 
accuracy  of  MR-based  DRRs  (characterized  by  the  MR-based  outlines/contour  lines). 
Points  3  and  4  are  used  for  the  L-R  position  setup,  points  1  and  5  or  2  and  6  are  used  for 
the  I-S  position  setup,  and  points  7  and  8  are  used  for  the  A-P  position  setup. 
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Abstract 

MR  imaging  based  treatment  planning  for  radiotherapy  of  prostate  cancer  is 
limited  due  to  MR  imaging  system  related  geometrical  distortions,  especially 
for  patients  with  large  body  sizes.  On  our  0.23  T  open  scanner  equipped  with  the 
gradient  distortion  correction  (GDC)  software,  the  residual  image  distortions 
after  the  GDC  were  <5  mm  within  the  central  36  cm  x  36  cm  area  for  a  standard 
48  cm  field  of  view  (FOV).  In  order  to  use  MR  imaging  alone  for  treatment 
planning  the  effect  of  residual  MR  distortions  on  external  patient  contour 
determination,  especially  for  the  peripheral  regions  outside  the  36  cm  x  36  cm 
area,  must  be  investigated  and  corrected.  In  this  work,  we  performed  phantom 
measurements  to  quantify  MR  system  related  residual  geometric  distortions 
after  the  GDC  and  the  effective  FOV.  Our  results  show  that  for  patients  with 
larger  lateral  dimensions  (>36  cm),  the  differences  in  patient  external  contours 
between  distortion-free  CT  images  and  GDC-corrected  MR  images  were 
1-2  cm  because  of  the  combination  of  greater  gradient  distortion  and  loss 
of  field  homogeneity  away  from  the  isocentre  and  the  uncertainties  in  patient 
setup  during  CT  and  MRI  scans.  The  measured  distortion  maps  were  used  to 
perform  point-by-point  corrections  for  patients  with  large  dimensions  inside 
the  effective  FOV.  Using  the  point-by-point  method,  the  geometrical  distortion 
after  the  GDC  were  reduced  to  <3  mm  for  external  contour  determination  and 
the  effective  FOV  was  expanded  from  36  cm  to  42  cm. 

(Some  figures  in  this  article  are  in  colour  only  in  the  electronic  version) 


1.  Introduction 

Computed  tomography  (CT)  has  played  a  very  important  role  in  radiation  treatment  planning 
in  modern  radiotherapy.  Based  on  the  information  provided  by  CT  data,  a  treatment  planning 
system  can  delineate  patients’  tumour  volumes,  external  contours  and  critical  anatomic 
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structures,  obtain  digitally  constructed  radiographs  (DDR)  for  patient  treatment  set-up  and 
perform  accurate  dose  calculations  for  treatment  planning.  Although  CT-based  treatment 
planning  has  been  considered  as  a  ‘gold  standard’  in  radiotherapy  for  years,  it  has  long  been 
recognized  as  having  poor  soft  tissue  contrast,  which  may  bring  uncertainties  into  treatment 
planning.  For  example,  Khoo  et  al  (1999)  found  that  MRI  provided  improved  definitions  of 
both  prostate  and  critical  structures  compared  to  CT.  Additionally,  Rasch  et  al  (1999)  found 
that  the  prostate  volume  was  40%  larger  on  CT  than  MR,  in  agreement  with  Roach  et  al  (1996), 
who  found  that  the  prostate  was  32%  larger  on  CT  than  on  MRI.  Debois  et  al  (1999)  showed 
that  improved  rectal  and  prostate  volume  delineation  from  MRI  could  lead  to  improvement  in 
both  target  coverage  and  rectal  sparing.  As  many  investigators  have  demonstrated  (Zelefsky 
et  al  1998,  Pollack  et  al  1999,  Pollack  et  al  2000,  Hanks  et  al  1998,  Hanks  1999),  dose 
escalation  with  three-dimensional  (3D)  conformal  and  recently  IMRT  increases  local  control 
while  reducing  complications  in  nearby  critical  structures.  However,  as  dose  levels  are 
increased,  the  precise  knowledge  of  tumour  volumes  and  the  accuracy  of  dose  calculation 
and  delivery  become  critical.  Magnetic  resonance  (MR)  imaging  has  been  shown  to  provide 
more  consistent  tumour  volume  delineation  in  treatment  planning  due  to  its  superior  soft  tissue 
contrast  than  CT  for  a  variety  of  sites. 

Although  CT  and  MR  fusion  has  been  widely  accepted  as  a  practical  approach  for  both 
accurate  delineation  (using  MR  data)  and  dose  calculation  (using  CT  data)  it  would  be  ideal 
if  MRI  could  replace  CT  entirely  for  treatment  planning,  retaining  the  superior  soft  tissue 
contrast  for  image  segmentation  and  at  the  same  time  eliminating  the  potential  errors  in  image 
fusion  due  to  (a)  variations  in  patient  set-ups  between  CT  and  MR  and  (b)  uncertainties  in  target 
localization  between  the  two  image  modalities.  Furthermore,  MRI-based  treatment  planning 
will  avoid  redundant  CT  imaging  sessions,  which  in  turn  will  avoid  unnecessary  radiation 
exposure  to  patients,  reduce  treatment  costs,  and  save  patient,  staff  and  machine  time.  There 
has  been  a  growing  interest  in  the  utilization  of  MRI-based  radiotherapy  treatment  planning  for 
the  past  decade,  and  progress  has  been  made  in  developing  an  effective  MRI-based  method  for 
treatment  planning  (Beavis  et  al  1998,  Guo  1998,  Mah  et  al  2002,  Chen  et  al  2002,  Michiels 
et  al  1994,  Chen  et  al  2003,  Lee  and  Bollet  2003). 

It  is  well  known  that  there  are  several  technical  challenges  for  MRI-based  treatment 
planning:  (1)  to  minimize  MR  image  geometrical  distortions,  (2)  to  incorporate  sequences 
that  reduce  MR  distortions,  (3)  to  develop  methods  that  can  correct  MR-related  distortions, 
(4)  to  correlate  MR  signal  intensity  with  the  densities  of  the  materials  imaged  and  (5)  to  obtain 
MRI-derived  DRR  for  patient  set-up. 

We  have  been  working  extensively  on  developing  MRI-based  treatment  planning  models 
for  the  last  few  years  (Chen  et  al  2003,  2004a,  2004b,  2004c).  It  has  been  generally  accepted 
that  there  are  no  clinically  significant  differences  in  dose  calculation  between  homogeneous 
and  heterogeneous  geometry  for  pelvis  and  it  has  been  common  practice  to  use  homogenous 
geometry.  We  have  implemented  MRI-based  prostate  treatment  planning  at  Fox  Chase  Cancer 
Center  (FCCC).  We  have  verified  the  dose  calculation  accuracy  with  unit-density  geometry 
built  for  MR  data.  The  agreement  was  within  2%  with  CT-based  dose  calculation  (Chen  et  al 
2004a).  An  effective  procedure  has  been  achieved  to  derive  MRI-based  DRRs  for  prostate 
cancer  patient  treatment  set-up  (Chen  et  al  2004b). 

Our  clinical  MR  unit  (0.23  T  open  MR  system.  Philips  Medical  Systems,  Cleveland, 
OH)  is  equipped  with  the  gradient  distortion  correction  (GDC)  software  (Mah  et  al  2002). 
Comparisons  of  patients’  MR  and  distortion-free  CT  images  showed  that  the  systematic 
geometrical  distortions  were  very  small  after  the  GDC  for  patients  with  36  cm  lateral 
dimensions  and  the  dosimetric  consequences  due  to  these  distortions  were  negligible  for 
treatment  planning.  However,  for  patients  with  larger  lateral  dimensions,  the  differences  in 
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Figure  1.  The  Philips  0.23  Tesla  open  MR  unit  (a)  and  the  FI 8  phantom  (b).  The  FI 8  phantom 
with  the  support  structure  is  placed  in  the  body /spine  coil.  The  support  structure  allows  the  marker 
plane  to  align  with  any  selected  imaging  planes. 


patient  external  contours  between  distortion-free  CT  images  and  GDC-corrected  MR  images 
were  significant  (1-2  cm)  because  of  the  combination  of  greater  gradient  distortion  and  loss 
of  field  homogeneity  away  from  the  isocentre  and  the  uncertainties  in  patient  set-up  during 
CT  and  MR1  scans.  The  uncertainties  on  the  external  contours  will  lead  to  errors  in  beam  path 
length  determination,  which  in  turn  may  result  in  dose  errors  of  clinical  significance.  In  order 
to  implement  MR-based  treatment  planning  clinically,  the  effective  FOV  of  our  MR  scanner 
must  be  investigated  and  expanded  to  include  patients  with  larger  lateral  dimensions. 

In  this  paper,  we  present  the  results  of  our  MR  image  distortion  investigation  on  a  0.23  T 
open  scanner  for  prostate  IMRT  treatment  planning.  A  point-by-point  mapping  technique  is 
developed  to  further  reduce  residual  gradient  distortions  after  GDC  post-processing.  We  will 
describe  the  technical  aspects  of  our  method  and  its  implementation  and  validation.  We  will 
show  that  together  with  set-up  uncertainty  quantification,  the  effective  FOV  of  the  scanner 
can  be  expanded  to  include  almost  all  the  prostate  patients,  and  the  uncertainties  of  patient 
external  contours  due  to  MR  image  distortions  can  be  determined  accurately. 

2.  Materials  and  method 

2.1.  The  MR  scanner 

The  MR  system  used  for  this  study  was  a  0.23  Tesla  open  scanner  operating  clinically  in  the 
department  of  Radiation  Oncology  at  Fox  Chase  Cancer  Center  (figure  1).  The  MR  scanner 
consists  of  two  approximately  1  m  diameter  poles.  The  separation  between  the  two  poles  is 
47  cm.  The  MR  scan  table  can  be  moved  in  orthogonal  planes  along  a  set  of  rails  mounted 
on  the  floor  and  on  an  orthogonal  set  of  rails  built  in  the  couch.  A  flat  table-top  made  of 
special  material,  which  is  stiff  and  light,  was  inserted  beneath  the  patient.  A  set  of  pads  made 
of  special  foam  was  used  to  adjust  the  table  height  according  to  the  patient  size.  The  three 
triangulation  lasers  (centre  and  laterals)  identical  to  those  used  on  the  linear  accelerators  were 
used  for  patient  positioning. 

2.2.  MR  geometrical  distortions  and  the  GDC  software 

In  MR  imaging,  each  point  in  the  imaging  space  is  associated  with  a  resonance  frequency, 
which  is  generated  by  the  linear  magnetic  field  gradient.  In  one  dimension,  say  the  x  direction, 
this  can  be  expressed  by 

m(x)  =  y(B0  +xGx) 


(1) 
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where  Bq  and  Gx  are  the  main  magnetic  field  intensity  and  the  magnetic  field  gradient  in 
the  x  direction  respectively,  and  y  is  the  gyromagnetic  ratio.  It  is  common  knowledge  that 
non-uniformity  of  Bq  and  nonlinearity  of  Gx  will  introduce  uncertainties  to  the  images  on  the 
spatial  localization  of  the  different  resonance  frequencies  because  of  the  way  MR  images  are 
reconstructed.  These  uncertainties  were  defined  as  MR  imaging  geometrical  distortions. 

MR  imaging  geometrical  distortions  can  be  divided  into  machine-induced  geometrical 
distortion  and  patient-induced  geometrical  distortion  (Michiels  et  al  1994).  The  machine- 
induced  MR  geometrical  distortion  includes  the  main  field  inhomogeneity,  gradient  field 
nonlinearity  and  eddy  currents  caused  by  the  gradient  switching.  As  has  been  demonstrated 
by  many  investigators  (Michiels  et  al  1994,  Chang  et  al  1990,  Sumanaweera  et  al  1994,  Schad 
et  al  1987),  the  distortions  induced  by  gradient  nonlinearity  and  main  magnetic  field  non¬ 
uniformity  are  stable  and  can  be  corrected  independently  of  other  machine  induced  distortions. 
The  patient-induced  geometrical  distortions  include  susceptibility  effects,  chemical  shift  and 
flow,  which  generally  cannot  be  corrected  ‘once  and  for  all’  because  of  its  patient  dependency. 
However,  due  to  the  high  receiver  bandwidth  (>100  Hz/pixel)  in  the  frequency  encoding 
direction  used  in  our  routine  scanning,  and  based  on  the  fact  that  there  are  no  patient-induced 
MR  image  artefacts  detected  with  0.2  T  low  magnetic  field  using  receiver  bandwidth  >  100  Hz/ 
pixel  in  the  frequency  direction  (Fransson  et  al  2001),  the  patient-induced  geometrical 
distortion  was  considered  negligible  for  our  MR  unit.  Therefore,  the  machine-induced 
geometrical  distortions  to  the  images  from  our  MR  scanner  can  be  detected,  studied  and 
corrected  accurately. 

A  Federal  Drug  Administration  (FDA)  approved  GDC  software  provided  by  the 
manufacturer  was  installed  on  our  scanner  for  MR  image  post-processing.  The  effectiveness 
of  the  post-processing  was  evaluated  using  a  standard  phantom  (FI 8  phantom)  provided 
by  the  manufacturer  (figure  1(b)).  The  F18  phantom  is  intended  for  evaluation  of  residual 
geometrical  distortions  in  the  images  post-processed  by  the  GDC  software.  The  phantom 
consists  of  a  rectangular  grid  of  spherical  markers  in  a  single  plane  separated  by  2.5  cm  in 
both  directions  (17  columns  and  15  rows).  The  phantom  is  3  cm  in  thickness  and  supplied 
with  a  support  structure  that  allows  the  F18  phantom  to  be  fixed  into  a  body/spine  coil.  The 
orientation  of  the  marker  plane  can  be  selected  to  align  with  any  of  the  fixed  imaging  planes. 
In  our  measurements,  the  F18  phantom  was  placed  in  the  xy  plane  of  the  coordinate  system 
(transverse  plane)  with  z  =  Zo,x  defined  as  horizontal  and  y  as  vertical. 

It  should  be  emphasized  that  MRI  geometrical  distortions  are  caused  by  non-uniformity 
of  the  magnetic  field  and  gradient  nonlinearity  in  the  three-dimensional  space.  The 
residual  distortions  after  the  GDC  are  expected  to  be  three-dimensional.  Since  MRI  does 
not  provide  electron  density  information  required  for  inhomogeneity  corrections  in  dose 
calculation  and,  for  some  sites  like  the  prostate,  these  corrections  are  not  clinically  significant 
(Photon  Treatment  Planning  Collaborative  Working  Group  1991),  we  have  used  homogeneous 
geometry  for  both  CT-based  and  MRI-based  treatment  planning.  The  homogeneous  geometry 
makes  the  distortion  correction  in  the  inferior-superior  direction  unnecessary  in  dose 
calculation.  The  only  corrections  needed  in  MRI-based  treatment  planning  are  those  in 
the  transverse  plane,  which  affect  the  patient  external  contours.  In  this  study,  only  distortions 
in  the  transverse  planes  were  considered. 

The  effective  FOV  after  the  GDC  was  established  36  cm  within  which  the  effect  of  residual 
MR  distortions  on  patient  external  contours  was  less  than  5  mm.  A  larger  effective  FOV  is 
required  for  MRI-based  treatment  planning  of  patients  with  larger  lateral  dimensions.  In  this 
work  we  have  quantified  the  residual  distortions  after  the  GDC  in  the  areas  outside  the  36  cm 
FOV  using  the  F18  phantom.  A  Synergy  body/spine  L  coil  was  selected  and  the  phantom 
(with  the  support  structure)  was  inserted  into  the  coil,  guided  by  the  laser  beam.  The  coil  was 
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placed  at  the  centre  of  the  0.23T  scanner  with  the  phantom  plane  in  the  sagittal  direction  (the 
x  and  y  plane  of  the  coordinate  system).  The  phantom  was  imaged  at  different  positions  with 
a  small  increment  ± A"  in  the  z-direction.  The  images  were  then  exported  in  the  DICOM 
format  and  transferred  to  our  development  environment  for  processing.  Since  each  marker’s 
geometrical  position  on  the  phantom  is  known  (also  can  be  calculated  accurately  knowing 
the  position  of  the  centre  of  the  phantom  on  the  image),  the  distortions  for  the  markers  at 
different  z  positions  can  be  quantified  by  comparison  with  the  known  geometrical  positions 
of  the  markers.  For  each  position  along  the  z-direction,  a  mapping  file,  which  records  the 
distortions  in  the  transverse  plane,  was  derived  based  on  the  measurements  of  the  markers, 
and  a  computer  program  was  developed  to  correct  the  distortions  point-by-point  for  any  MR 
scans  using  the  measured  distortion  maps  stored  in  the  files. 

Since  the  FI 8  only  has  17  markers  in  the  horizontal  direction  and  15  in  the  vertical 
direction,  the  distortions  outside  the  40  x  35  cm2  central  area  were  measured  by  shifting  the 
phantom  left  and  right  a  few  centimetres  to  determine  the  effective  FOV. 

2.3.  The  point-by-point  distortion  correction  method 

Based  on  the  distorted  images  of  the  grid  points  at  different  z  positions,  correction  maps  can 
be  derived  by  comparisons  of  the  distorted  images  with  the  physical  pattern  of  the  grid  on 
the  FI 8  phantom.  These  maps  can  then  be  applied  to  rectify  patient  images  acquired  at  the 
corresponding  positions,  using  identical  imaging  sequences  and  parameter  settings. 

To  quantify  the  distortion,  i.e.  the  differences  between  the  corresponding  marker  points  of 
the  F 1 8  phantom  and  the  image  acquired  from  the  scanner,  a  Cartesian  coordinate  is  introduced, 
and  the  origin  (0,  0)  of  the  coordinate  is  set  at  the  bottom  left  comer  of  the  images.  In  the 
Cartesian  coordinate,  we  define  the  distributions  of  the  grid  points  as  r? (i,  j )  and  rs„(i,  j), 
where  r? (/',  j)  is  the  spatial  distribution  of  the  points  on  the  FI 8  phantom  and  ri(i,  j)  is  the 
spatial  distribution  of  the  centres  of  the  markers  on  the  image  scanned,  respectively.  The 
distortion  maps  can  then  be  defined  as 

MAp  J)  =  rP(i,  j)  -  r'z(i,  j)  (2) 

where  (=1,2,...  ,17  and  j  —  1,2,...  ,15. 

Once  the  maps  of  the  geometrical  distortions  of  the  markers  are  established,  they  can  be 
used  to  correct  the  residual  distortions  of  MR  images. 

Figure  2  demonstrates  the  principles  of  mapping  the  pixels  within  the  area  defined  by  four 
marker  points  on  the  FI 8  phantom  and  the  corresponding  area  on  the  distorted  images.  We 
assumed  that  the  mapping  area  on  the  distorted  image  is  quadrilateral,  which  may  introduce 
uncertainties  in  the  distortion  corrections  since  a  quadrangle  may  not  have  a  quadrilateral 
image  due  to  gradient  distortion  and  non-uniformity  of  the  main  magnetic  field.  Flowever,  the 
quadrilateral  approach  is  a  reasonable  one  since  it  will  provide  the  required  accuracy  for  our 
application. 

Based  on  the  quadrilateral  approach,  we  can  perform  the  point-by-point  correction  to 
the  distorted  images.  First,  those  pixels  inside  the  area  defined  by  any  four  markers  on  the 
phantom  image  are  correlated  with  these  on  the  distorted  image  by  a  linear  mapping.  This 
is  done  by  equally  dividing  the  distance  between  any  two  points  on  both  images  and  then 
correlating  the  areas  defined  by  the  connecting  lines. 

The  area  (distorted  pixels)  on  the  distorted  image  may  be  larger  or  smaller  than  the 
area  (pixels)  on  the  F18  phantom.  Figure  3  shows  a  distorted  pixel  and  its  relationship  with 
pixels  on  the  distorted  MR  image.  Assume  that  a  pixel  ( m ,  n )  on  the  F18  phantom  has  its 
corresponding  distorted  area  on  the  MR  image  that  is  larger  than  the  MR  pixel  size  as  shown 
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Figure  2.  A  diagram  showing  the  mapping  of  the  pixels  between  the  distorted  image  and  the  FI 8 
phantom. 
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Figure  3.  A  diagram  to  illustrate  the  pixel  intensity  calculation. 


in  figure  3.  The  area  covers  not  only  the  pixel  ( i ,  j )  but  also  part  of  pixels  (i  —  1,  j),  (i,  j  —  1), 
O',  j  +  1),  0  +  1 ,  j  +  1)  and  0  +  1,  j)-  The  overlapping  areas  here  are  represented  by  A(k,  l) 
with  k  =  i  —  1,  i,  i  +  1  and  l  =  j  —  I ,  /.  /  +  I .  Then,  the  corrected  intensity  value  for  pixel 
(m,  n)  can  be  calculated  as 

/(m,«)  =  ^A(M)'/'(U)  (3) 

k,l 

where  k  —  i  —  1,  i,  i  +  1  and  l  =  j  —  1,  j,  j  +  1,  S  represents  the  value  of  the  area  (m,  n) 
and  I'(k,  l )  is  the  image  intensity  for  pixel  ( k ,  Z)  on  the  distorted  MR  image.  If  the  distorted 
area  of  a  pixel  is  smaller  than  the  size  of  a  pixel,  the  image  intensity  can  be  approximated  by 
choosing  I(m,  n )  =  /'(£,  /). 

2.4.  Evaluation  of  the  point-by-point  distortion  correction  method 

Since  there  is  no  patient-induced  distortion  for  the  0.23  T  MR  unit  with  the  sequences  used  in 
our  study,  the  distortion  maps  derived  from  the  phantom  measurement  can  be  used  to  correct 
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Figure  4.  Actual  viewable  areas  of  the  scanner  and  a  distortion  map.  (a)  The  outer  contour 
represents  the  viewable  area  of  the  transverse  plane  at  the  isocentre  of  the  scanner  and  the  imier 
contour  the  viewable  area  at  the  transverse  planes  10  cm  away  from  the  isocentre  plane  along 
the  I-S  direction.  The  viewable  areas  for  any  transverse  planes  within  10  cm  of  the  isocentre 
will  be  between  the  two  contours,  (b)  A  distortion  map  measured  based  on  the  FI 8  phantom 
measurement.  Without  distortion,  any  adjacent  marker  points  laterally  and  vertically  would  have 
a  2.5  cm  separation. 


the  residual  MR  distortions.  However,  detailed  testing  on  real  patients  must  be  performed 
before  its  clinical  implementation.  In  this  work,  distortion-free  CT  images  were  used  to 
compare  with  MRI  images  of  prostate  patients  before  and  after  the  GDC  corrections. 

It  was  understood  that  the  differences  between  the  CT  and  MRI  external  contours  for  the 
same  patients  resulted  not  only  from  MR  image  distortions  but  also  from  set-up  uncertainties 
between  CT  and  MR  scans,  which  had  to  be  minimized  or  excluded  in  order  to  evaluate  the 
accuracy  of  the  distortion  correction.  As  previous  investigators  reported  (Mah  et  al  2002, 
Chen  et  al  2005),  the  set-up  uncertainties  in  CT  and  MR  could  be  about  0.5  cm  respectively. 
The  combined  set-up  uncertainties  in  our  CT  and  MRI  comparisons  were  therefore  expected 
to  be  ~0.7  cm  (one  standard  deviation).  Clinically,  we  have  observed  set-up  uncertainties 
between  CT  and  MR  scans  larger  than  1  cm  for  some  patients  with  large  lateral  dimensions. 

In  order  to  evaluate  our  point-by-point  distortion  correction  method,  static,  porcine  tissue 
samples  were  used  in  this  study  to  eliminate  the  set-up  uncertainties.  We  imaged  frozen, 
porcine  leg  quarters  on  our  CT  and  MRI  scanners  and  arranged  a  pair  of  these  to  mimic  human 
pelvis,  using  a  special  cradle  to  hold  their  positions  firmly.  All  the  MR  images  were  scanned 
with  our  routine  clinical  three-dimensinal  fast  spin  echo  sequences  (3DFSE,  256  x  256, 
1.95  mm  pixel,  3  mm  slice  thickness,  TR  =  3000  ms,  TE  =  140  ms,  BW  readout  gradient 
>100  Hz/pixel),  and  a  standard  FOV  of  48  cm.  The  external  contours  from  the  CT  images 
were  exported  to  the  MR  images  after  fusion  for  direct  comparison. 


3.  Results  and  discussion 

Figure  4(a)  shows  the  actual  viewable  areas  of  our  0.23  T  MR  unit  for  a  standard  48  cm 
FOV  for  prostate  imaging  that  indicate  the  geometrical  limitations  of  the  scanner  at  different 
transverse  planes.  The  marker  points  outside  the  viewable  areas  were  outside  the  256  by  256 
matrix  due  to  the  MR  geometric  distortions  and  were  lost  in  the  initial  raw  MR  images.  The 
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(a)  (b)  (c) 

Figure  5.  MR  images  of  the  FI 8  phantom  at  the  isocentric  plane:  (a)  before  the  GDC,  (b)  after 
the  GDC,  and  (c)  further  corrected  by  the  point-by-point  distortion  correction. 


Figure  6.  CT  and  MR  images  of  a  prostate  patient:  (a)  a  distortion-free  CT  image,  (b)  the  same 
patient  imaged  on  the  MR  scanner  without  any  distortion  corrections,  (c)  only  with  the  GDC  post¬ 
processing,  and  (d)  with  both  the  GDC  and  our  point-by-point  correction.  The  external  contours 
on  the  MR  images  were  exported  from  the  CT  image  after  image  fusion. 


GDC  software  reduced  the  MR  geometric  distortions  significantly  but  could  not  recover  the 
missing  marker  points.  The  outer  contour  in  figure  4(a)  represents  the  actual  viewable  area 
of  the  scanner  at  the  isocentre  transverse  plane  and  the  inner  contour  represents  the  actual 
viewable  area  which  is  10  cm  away  from  the  isocentre.  Any  transverse  planes  within  10  cm 
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Figure  7.  CT  and  MR  images  of  frozen  porcine  samples:  (a)  a  distortion-free  CT  image  at  the 
isocentric  plane,  (b)  the  same  isocentric  slice  imaged  on  the  MR  scanner  without  any  distortion 
corrections,  (c)  only  with  the  GDC  post-processing  and  (d)  with  both  the  GDC  and  our  point-by¬ 
point  correction.  The  external  contours  on  the  MR  images  were  exported  from  the  CT  image  after 
image  fusion. 


of  the  isocentre  will  have  a  viewable  area  between  the  two  contours.  The  actual  viewable  area 
decreases  with  increasing  distances  from  the  isocentre,  which  is  expected  because  of  the  loss 
of  uniformity  of  the  magnetic  field  at  large  distances.  The  distortion  maps  were  established  for 
every  transverse  plane  along  the  I-S  direction  with  a  0.5  cm  interval  and  the  maps  were  used  to 
correct  the  distortions  in  our  point-by-point  distortion  correction  method.  Figure  4(b)  shows 
the  distortion  map  of  a  transverse  plane  at  the  location  of  the  scanner’s  centre.  Compared 
with  the  physical  pattern  of  the  FI 8  phantom,  the  residual  distortions  after  the  GDC  can  be 
observed  and  the  maximum  distortion  is  about  1  cm. 

Figure  5  shows  MR  images  of  the  FI  8  phantom  with  and  without  the  GDC  post-processing 
and  after  the  point-by-point  corrections.  These  images  were  taken  at  the  centre  of  the  scanner 
with  a  standard  FOV  of  48  cm  for  prostate  imaging.  We  then  established  a  new  effective  FOV 
based  on  the  viewable  area  for  this  transverse  plane,  which  showed  clinically  insignificant 
distortions  after  the  point-by -point  correction.  We  established  the  effective  FOVs  for  all  the 
transverse  planes  in  the  same  way.  For  later  routine  clinical  applications,  only  the  image 
points  within  the  effective  FOVs  were  processed.  The  maximum  distortions  were  about 
1.0  cm  before  the  point-by-point  correction  and  were  reduced  to  <0.3  cm  after  the  correction. 
The  improvement  on  the  marker  positions  is  shown  clearly. 

A  few  patients  were  selected  and  processed  to  test  our  method.  We  found  that 
the  distortions  were  reduced  dramatically  after  applying  the  GDC  and  our  point-by-point 
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corrections.  Figure  6  shows  a  distortion-free  CT  image  with  an  external  contour  drawn  by  a 
physician  and  two  MR  images  with  and  without  the  GDC  and  a  MR  image  with  both  the  GDC 
and  the  point-by-point  distortion  correction.  The  maximum  residual  error  after  our  point-by- 
point  distortion  correction  is  expected  to  be  ~0.3  cm  based  on  the  investigation  results  using 
the  F18  phantom  (up  to  2  cm  distortions  were  corrected  compared  with  the  initial  image). 
From  figure  6(d),  the  differences  between  the  CT  external  contour  and  the  distortion  corrected 
image  are  still  seen.  It  was  believed  that  the  major  reason  for  these  differences  was  the  set-up 
uncertainty  between  the  CT  and  MR  scans. 

To  demonstrate  this,  we  performed  paired  CT  and  MR  scans  for  large,  static  porcine 
samples  immobilized  in  a  specially  made  cradle.  Figure  7(a)  shows  a  distortion-free  CT 
image  of  the  porcine  samples,  which  were  42  cm  in  lateral  dimensions.  The  external  contour 
of  the  CT  image  was  transferred  to  the  MR  images  for  comparison.  The  maximum  distortion 
measured  on  the  initial  MR  image  was  2.1  cm  (figure  7(b))  and  significant  residual  distortions 
were  still  observable  after  the  DGC  post-processing  (figure  7(c)).  Measurements  on  the  final 
MR  image  showed  that  the  MR  geometrical  distortions  were  reduced  to  <0.3  cm  after  our 
point-by-point  distortion  correction  (figure  7(d)). 

4.  Conclusions 

In  this  paper,  we  have  described  our  work  on  MR  image  distortion  correction  to  further 
improve  the  accuracy  of  dose  calculation  for  MR-based  treatment  planning  for  prostate  cancer. 
Our  goal  was  to  develop  practical  methods  for  the  clinical  implementation  of  MRI-based 
treatment  planning.  Our  studies  showed  that,  with  our  routine  clinical  three-dimensinal  fast 
spin  echo  sequences  (3DFSE,  256  x  256,  1.855  mm  pixel,  TR  =  140  ms,  TE  =  3000  ms, 
BW  readout  gradient  >100  Hz/pixel),  there  were  no  patient-induced  distortions.  Therefore, 
the  residual  machine  specific  geometrical  distortions  after  the  GDC  could  be  quantified  by 
phantom  measurements  and  further  reduced  by  our  point-by-point  correction  technique.  The 
effective  FOVs  of  the  scanner  were  established  based  on  the  actual  viewable  areas  with 
adequate  geometric  distortion  corrections  (ensuring  <5  mm  distortion  error).  The  effective 
FOV  for  prostate  imaging  using  a  standard  FOV  of  48  cm  has  been  expanded  from  36  cm  using 
the  existing  GDC  software  to  42  cm  using  the  point-by-point  distortion  correction  technique 
developed  in  this  work.  Our  results  indicated  that,  with  the  distortion  maps  established  in  this 
work,  we  could  correct  MR  geometrical  distortions  for  patients  of  lateral  dimensions  up  to 
42  cm.  Significant  improvement  in  dose  calculation  has  been  achieved  based  on  a  1-2  cm 
improvement  in  patient  external  contour  determination. 
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MR-based  treatment  planning  for  prostate  IMRT  has  been  successfully  developed  and  implemented  in  our  department.  In  a 
previous  work,  we  have  performed  dose  verification  on  a  commercial  treatment  planning  system  (TPS)  for  the  MR-based 
treatment  planning.  However,  due  to  the  limitation  of  our  low  field  MR  scanner  (2  mm  achievable  resolution)  and  the 
inflexibility  of  image  processing  on  the  TPS,  uncertainties  were  existed  in  the  comparisons  of  the  CT  and  the  MR-based  plans. 
Our  previous  work  also  showed  the  discrepancies  in  dose  calculations  in  IMRT  planning  between  using  Monte  Carlo  method 
and  the  TPS  with  non-coplanar  beam  arrangements  that  a  large  amount  of  bones  in  pelvic  area  were  irradiated.  In  this  work,  we 
used  the  Monte  Carlo  to  further  verify  the  dosimetric  accuracy  and  consistency  for  the  MR-based  IMRT  treatment  planning 
and  to  study  the  inhomogeneity  effect  by  using  bulk  density  assignments  to  bony  structures  on  MR  images.  Fifteen  prostate 
patients  were  studied  using  homogeneous  geometry  based  on  CT  and  MRI,  and  heterogeneous  geometry  built  based  on  CT 
numbers  or  different  bulk  densities  for  MR  bony  structures.  For  IMRT  treatments  with  coplanar  beam  arrangements,  the  mean 
dose  values  of  the  CTV  for  homogeneous  geometry  based  on  CT  were  about  2%  higher  than  those  for  heterogeneous  geometry 
based  on  CT.  The  difference  in  the  mean  CTV  dose  between  homogenous  MRI  and  heterogeneous  CT  geometries  was  about 
3%  and  less.  The  DVH  curves  also  agreed  to  within  5%  in  dose  or  volume  among  these  plans.  For  treatments  with  larger 
discrepancies  found  in  the  CTV  dose  between  homogeneous  and  heterogeneous  geometry  due  to  the  attenuation  of  the  oblique 
beams  going  through  the  femurs,  bulk-density  assignment  to  femurs  and  femoral  heads  were  used  and  a  1.8g/cm3  optimal 
bulk-density  was  found  to  reduced  the  discrepancies  from  10%  to  less  than  3%  in  MR-based  planning. 

I.  Introduction 

Since  the  past  a  couple  of  years,  a  sophisticated  technique  and  a  practical  procedure  have  been  developed  and  implemented 
clinically  in  our  department  for  MR-based  treatment  planning  for  prostate  cancer.11 !.  We  have  developed  an  efficient  method  to 
quantify  and  correct  the  machine  induced  geometrical  distortions  to  less  than  0.3  cm.121  We  have  performed  dose  verification 
for  MR-based  treatment  planning  on  a  commercial  treatment  planning  system  (TPS)  and  good  agreements  were  achieved 
between  CT  and  MR-based  dose  calculations. !  1 !  However,  due  to  the  limitation  of  our  low  field  MR  scanner  (2  mm  achievable 
resolution)  and  the  inflexibility  of  image  processing  on  the  TPS,  some  discrepancies  were  found  in  our  previous  work  during 
the  comparisons  of  the  CT  and  the  MR-based  plans,  especially  for  the  doses  to  the  critical  structures,  such  as  doses  to  rectum 
and  bladder.  Discrepancies  in  dose  calculations  of  IMRT  planning  were  also  existed  between  using  Monte  Carlo  method  and 
the  TPS  with  non-coplanar  beam  arrangements  that  a  large  amount  of  bones  in  pelvic  area  were  irradiated.'3.  In  this  work,  we 
use  the  Monte  Carlo  method  to  perform  dose  calculation  to:  1)  further  verify  the  accuracy  and  consistency  of  the  dose 
calculation  in  the  MR-based  treatment  planning;  2)  investigate  the  inhomogeneity  effect  in  the  CT-based  calculation  and  to 
compare  the  CT  results  with  bulk-density  assignment  in  MR-based  calculations. 

II.  Materials  and  methods 

A.  IMRT  planning 

The  images  for  all  the  patients  in  this  work  were  acquired  on  a  Picker  AcQSim  CT  scanner  and  a  Philips  MR  0.23  T  open 
scanner  (Philips  Medical  System,  Cleveland,  OH).  Target  and  critical  structures  were  delineated  by  the  physicians  based  on  the 
fused  CT  and  MR  images  using  the  AcQSim  VoxelQ  segmentation  tools.  The  IMRT  plans  were  generated  on  the  Corvus 
inverse-planning  system  (Version  5)  based  on  the  CT  data  (Nomos  Corp.  Sewickley,  PA).  Once  the  plan  was  approved  by  the 
physician,  it  was  moved  to  our  Monte  Carlo  system  with  the  patient’s  CT  and  contour  data,  treatment  isocenter,  leaf¬ 
sequencing  files  and  beam  parameters  for  dose  calculation.  On  the  Monte  Carlo  treatment  planning  system,  images  were 
further  processed  so  that  the  FOV,  voxel  size  and  volumes  of  RT  structures  were  kept  the  same  for  both  CT-based  and  MR- 
based  dose  calculations. 

B.  Monte  Carlo  dose  calculation 

The  Monte  Carlo  code  used  in  this  work  for  dose  calculation  was  MCSIM,  which  is  an  EGS4/PRESTA  user  code  developed  at 
Fox  chase  Cancer  Center  (FCCC).'4'5'  The  beam  information  was  represented  using  a  source  model,  which  was  built  based  on 
measured  beam  data'8’7’6*,  and  validated  for  Monte  Carlo  dose  calculation  for  photon  beams  from  our  Siemens  accelerators. 
During  the  calculation,  the  multi-leaf  collimator  leakage  effect  was  taken  into  account  when  intensity  maps  were  reconstructed 
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from  a  plan.  The  accuracy  of  the  dose  calculation  is  better  than  2%  compared  with  measured  data.'9'101  For  each  patient,  an 
RTP  file  from  the  Corvus  treatment  planning  system  that  includes  patient  setup  parameters  and  beam  and  leaf-sequence 
information  was  used  for  the  Monte  Carlo  simulation.  The  electron  and  photon  energy  cutoffs,  ECUT  and  PCUT,  for  particle 
transport  in  the  Monte  Carlo  simulations  were  set  to  700  keV  and  10  keV,  respectively.  The  energy  thresholds  for  5  ray 
production  (AE)  and  for  Bremsstrahlung  production  (AP)  were  also  set  to  700  keV  and  10  keV,  respectively.  The  maximum 
fractional  energy  loss  per  electron  step  (ESTEP)  was  set  to  0.04  and  the  default  parameters  were  used  for  the  PRESTA 
algorithm.  Seventy  million  particle  histories  were  used  to  achieve  less  than  0.5%  statistical  uncertainties  to  the  target  dose  for 
all  the  IMRT  plans.  Each  photon  was  split  20  times  to  improve  the  simulation  efficiency  using  the  photon-splitting  technique 
implemented  in  MCSIM. 

C.  Investigation  of  inhomogeneity  effect  for  MR-based  calculation 

CT-based  IMRT  Monte  Carlo  dose  calculations  with  and  without  heterogeneity  correction  were  performed  for  each  selected 
plan  in  order  to  investigate  heterogeneity  effect  caused  by  different  beam  angle  arrangements.  Based  on  the  results,  MR-based 
IMRT  dose  calculations  were  then  performed  to  the  plan  using  either  uniform  density  geometry  or  uniform  density  geometry 
with  bulk  electron  density  assigned  to  bony  structures.  For  the  plans  with  insignificant  inhomogeneity  effect,  uniform 
geometries  with  water  density  were  used  in  the  MR-based  dose  calculation.  For  the  plans  that  bony  structure  constitutes  a  large 
part  of  volume  irradiated,  uniform  density  geometry  with  bulk  electron  density  assigned  to  bony  structures  was  used  in  the 
MR-based  dose  calculation. 

D.  Dose  verification  of  MR-based  plan 

The  differences  in  isodose  distributions  and  dose -volume  histograms  (DVF1)  were  used  to  compare  the  results  of  the  CT-based 
and  the  MR-based  dose  calculations  for  each  IMRT  plan.  The  clinical  quantities  such  as  target  volume  (CTV),  the  mean  dose, 
maximum  and  minimum  dose  received  by  the  CTV  and  the  critical  structures  were  chosen  for  the  comparison.  The  maximum 
dose  was  defined  as  the  highest  dose  received  by  1%  of  the  target  volume  and  the  minimum  dose  was  defined  as  the  lowest 
dose  received  by  99%  of  the  target  volume,  respectively.  The  paired  CT  and  MR  data  for  any  patients  in  this  work  were  pre- 
processed  to  have  the  same  pixel  resolution.  A  special  computer  code  was  developed  to  convert  the  patient  CT  and  MR  image 
data  from  the  DICOM  format  to  geometries  specially  formatted  for  the  MCSIM  code. 

III.  Results  and  discussions 

With  proper  geometrical  distortion  corrections  for  the  MR  images  and  the  use  of  uniform  water  equivalent  density  geometry  in 
the  dose  calculations,  it  was  expected  that  the  IMRT  plans  based  on  CT  and  MR  data  would  agree  well.  Figure  1  shows  the 
DVH  curves  from  the  CT  and  MR-based  dose  calculations  for  one  of  the  10  IMRT  plans  investigated.  This  IMRT  plan  used  10 
MV  photon  beams  with  9  gantry  angles  on  a  Siemens  Primus  accelerator  with  39  fractions  to  deliver  78  Gy.  The  difference  in 
the  mean  dose  of  the  CTV  was  0.3%  between  the  CT-based  and  MR-based  calculations.  The  differences  in  the  maximum  dose 
and  in  the  minimum  dose  of  the  CTV  were  0.8%  and  0.2%  respectively.  Figure  2  shows  the  isodose  distributions  from  the  CT- 
and  MR-based  calculations  for  this  patient. 

Table  1  lists  the  differences  of  the  clinical  quantities  in  the  dose  calculations  using  MR  data  and  CT  data  with  homogeneous 
water  equivalent  geometries  for  10  IMRT  prostate  plans  with  co-planar  beam  arrangements,  as  well  as  the  differences  of  the 
clinical  quantities  between  MR-based  dose  calculations  with  uniform  water  equivalent  geometry  and  CT-based  dose 
calculations  with  heterogeneity  corrections.  The  data  presented  in  Table  1  were  calculated  as  the  percentage  differences 
between  CT-based  and  MR-based  results,  i.e.,  (D(T  -DMR)/DMR  ><100%,  where  Dcr  is  the  dose  from  CT-based  calculation  and 
Dmr  is  the  dose  predicted  by  MR-based  calculation.  ADmin  is  the  difference  in  the  minimum  dose  of  the  target,  ADmean  is  the 
difference  in  the  averaged  target  dose,  ADmax  is  the  difference  in  the  maximum  target  dose,  and  ADS  and  AD95  are  the 
differences  of  the  dose  received  by  5%  and  95%  of  the  target  volume,  respectively.  Table  1  also  shows,  with  the  co-planar 
beam  arrangements,  the  averaged  differences  using  MR  data  with  uniform  water  equivalent  geometries  and  CT  data  with 
heterogeneity  corrections  for  all  the  clinical  quantities  used  in  this  comparison  were  less  than  2%  ,  and  the  maximum 
differences  were  also  less  than  4%  and  the  averaged  differences  over  the  10  IMRT  plans  were  less  then  1.6%  for  all  the 
quantities  in  the  comparison,  indicating  that  the  uniform  geometry  was  a  good  approach  with  our  commonly  used  beam 
arrangements. 

For  some  clinical  cases  non-coplanar  beam  arrangements  where  a  large  amount  of  pelvic  bones  were  irradiated  in  the 
treatment,  heterogeneity  correction  became  necessary  in  planning.  More  than  8%  differences  in  dose  calculations  were 
observed  in  the  comparisons  of  DVFIs  for  some  CT-based  IMRT  plans  with  and  without  heterogeneity  correction.  To  utilize 
MR-based  planning  for  the  treatments  with  large  amount  of  pelvic  bones  irradiated,  inhomogeneity  corrections  must  be  taking 
into  account  and  bulk-density  assignment  to  bony  structures  was  proposed  since  there  is  no  simple  point-to-point  correlation 
between  MR  signal  intensities  and  electron  densities  of  the  materials  imaged'111.  Various  bulk-density  assignments  to  the 
femurs  and  femoral  heads  were  studied  in  this  work.  Figure  3  shows  the  changes  in  the  DVF1  with  different  bulk  density 
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assigned  to  the  femurs  (between  1.5  and  2.2  g/cnr’)  for  two  IMRT  plans  with  large  amount  of  pelvic  bones  irradiated.  Table  2 
shows  the  target  doses  based  on  CT  data  with  heterogeneity  correction  and  based  on  MR  data  with  and  without  bulk  densities 
assigned  to  the  femurs  and  femoral  heads  for  a  7-field  IMRT  plan  with  7-field  non-coplanar  beams.  The  averaged  target  doses 
for  each  beam  as  calculated  using  CT  data  with  heterogeneous  correction  and  using  MR  data  with  bulk-density  assignment 
agreed  well  compared  with  CT  uniform  water  equivalent  geometry.  For  those  beams  with  gantry/table  angles  as  275/340,  85/20 
and  85/0,  the  differences  in  the  average  target  doses  were  decreased  from  about  10%  with  uniform  water  equivalent  geometry 
to  about  3%  or  less  after  using  1.8g/cm3  bulk  density  assignment  to  femurs  and  femoral  heads.  The  changes  in  DVHs  using  CT 
data  with  1.8g/cm3  bulk  density  for  the  femurs  and  femoral  heads  also  confirmed  our  findings. 


Figure  1.  Comparison  of  the  DVHs  for  a  co-planar 
plan  calculated  based  on  CT  data  with  and  without 
heterogeneity  correction  and  based  on  MR  data  in 
uniform  geometry  (1  g/cm3). 


Figure  2.  Isodose  comparison  between  CT-based  dose  calculation 
without  heterogeneity  correction  and  MR-based  dose  calculation  in 
uniform  geometry  (1  g/cm3). 


The  values  of  1.8g/cm3  for  optimal  bulk-density  assignment  for  the  femurs  and  femoral  heads  were  obtained  based  on  the 
results  of  our  calculations  only.  We  have  used  different  bulk-density  assignments  in  the  dose  calculations  for  IMRT  plans,  4 
beam  conformal  plans  and  single  beams  with  different  incident  angles  and  compared  the  results  with  those  predicted  by  CT 
data  with  heterogeneity  corrections;  1.8g/cm3  gave  the  best  fit. 

For  critical  organs  in  prostate  treatment,  such  as  rectum  and  bladder,  the  dose  differences  between  CT-  and  MR-based 
calculations  were  negligible  compared  with  the  differences  in  the  target  dose  for  all  the  plans  studied  (see  figures  1  and  3).  The 
results  '  '  "  "  calcif  '  . . 


Figure  3.  Comparison  of  DVHs  for  non-coplanar  plans  between  CT-based  dose  calculation  with  heterogeneity  correction  and  MR-based  dose 
calculation.  Different  bulk  densities  were  assigned  to  the  femurs  in  MR-based  dose  calculations,  a)  a  7  beam  plan  with  table  angles  at  25°,  20°,  0°, 
0°,  0°,  340°,  and  330°;  and  b)  a  9  beam  plan  with  table  angle  at  25°,  20°,  0°,  0°,  0°,  0°,  0°,  340°  and  335°. 
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Table  1  Target  dose  differences  between  CT-based  and  MR-based  treatment  plans.  The  values  at  the  left  in  the  brackets  are  the  differences  in 
target  doses  calculated  based  on  MR  and  CT  with  water  equivalent  uniform  geometry;  the  values  at  the  right  in  the  brackets  are  the  differences 
in  target  doses  calculated  based  on  MR  with  water  equivalent  uniform  geometry  and  CT  with  heterogeneity  correction  for  10  IMRT  plans  of 
prostate  cancer  patients 
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Table  2  Differences  in  the  target  dose  with  and  without  heterogeneity  correction  and  with  bulk  density  assignment  for  a  7-field,  non-coplanar 
IMRT  plan.  A\%  represent  the  percentage  differences  with  and  without  heterogeneity  correction;  A2V0  represent  the  percentage  differences  with 
heterogeneity  correction  and  1.8  g/cm  bulk  density  assigned  to  femurs  and  femoral  heads. 
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IV.  Conclusion 

We  found  that  the  differences  in  the  results  between  CT  data  (with  heterogeneity  correction)  and  MR  data  (with  uniform  water 
equivalent  geometry)  were  about  3%  or  less  and  less  than  2%  in  the  mean  values  for  the  10  plans  with  beams  arranged  in  one 
axial  plane  (see  table  1).  For  MR-based  calculations  (with  homogeneous  geometry),  our  results  demonstrated  that  the 
differences  between  MR-based  calculations  and  CT-based  calculations  (without  heterogeneity  correction)  were  less  than  2% 
for  the  individual  patients  and  about  1%  in  the  mean  values,  which  proved  that  MR-based  IMRT  plans  can  be  used  to  replace 
CT-based  planning  clinically  (see  table  1).  The  1%  -  2%  differences  in  dose  calculations  were  mainly  caused  by  the  setup 
uncertainties  of  the  two  imaging  modalities  if  geometrical  distortions  on  the  MR  images  were  corrected  to  less  than  3  mm. 

For  treatments  in  which  relatively  large  amount  of  bones  are  irradiated,  MR-based  treatment  planning  with  homogeneous 
geometry  would  not  be  appropriate  because  of  the  excessive  attenuation  of  the  photon  beams  passing  through  bony  structures. 
Flowever,  by  assigning  bulk  densities  to  the  bony  structures  especially  for  the  femurs  and  femoral  heads,  the  dose  differences 
could  be  reduced  to  less  than  3%  (see  table  2).  The  bulk  densities  assigned  to  the  femurs  that  gave  the  best  fits  to  CT-based 
calculations  with  heterogeneity  correction  were  1.8g/cm3  in  our  simulations. 

We  have  established  a  practical  procedure  for  MR-based  treatment  planning,  in  which  MR  geometrical  distortions  are  properly 
corrected,  a  useful  DRR  technique'1’  is  developed/implemented  and  homogeneous  geometry  is  used  in  the  dose  calculation. 
Our  results  show  that  it  is  accurate  and  consistent  in  dose  calculation  for  prostate  cancer  treatment.  In  our  department,  one  third 
of  radiotherapy  treatments  are  for  prostate  cancer  and  homogeneous  geometry  has  been  used  in  prostate  planning  dose 
calculation.  By  using  MR-based  treatment  planning  the  savings  in  staff  time,  patient  time  and  treatment  costs  would  be 
significant  and  the  additional  radiation  exposure  from  CT  scans  can  be  avoided.  It  is  concluded  that  MR-based  radiotherapy 
treatment  planning  can  be  widely  utilized  clinically  for  prostate  cancer  treatment.  Further  investigations  are  carried  out  for 
other  sites  where  heterogeneity  effect  is  not  significant. 
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